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ABSTRACT 

We present new Chandra ACIS-S3 observations of Cassiopeia A which, when combined with earher 
ACIS-S3 observations, show evidence for a steady ^ 1.5-2% yr~^ decUne in the 4.2-6.0 keV X-ray 
emission between the years 2000 and 2010. The computed flux from exposure corrected images over the 
entire remnant showed a 17% decline over the entire remnant and a sUghtly larger (21%) dechne from 
regions along the remnant's western hmb. Spectral fits of the 4.2-6.0 keV emission across the entire 
remnant, forward shock filaments, and interior filaments indicate the remnant's nonthermal spectral 
powerlaw index has steepened by about 10%, with interior filaments having steeper powerlaw indices. 
Since TeV electrons, which give rise to the observed X-ray synchrotron emission, are associated with 
the exponential cutoff portion of the electron distribution function, we have related our results to a 
change in the cutoff energy and conclude that the observed decline and steepening of the nonthermal 
X-ray emission is consistent with a deceleration of the remnant's ~5000 km s~^ forward shock of « 
10-40 km s-i yr-i. 

Subject headings: acceleration of particles — ISM: individual objects (Cassiopeia A) — radiation 
mechanisms: non-thermal 



1. INTRODUCTION 

Supernova remnants (SNRs) have long been considered 
to be the primary source of Galactic cosmic-rays (CRs) 
below the knee of the cosmic-ray spectrum, ~ 10^^ eV. 
TeV 7-ray observations of SNRs such as RX J1713. 7-3946 
and RX J0 852. 0-4622 provide evide nce for the accelera- 
tion of ions (EE aronian et al.ll2007al lbl). However the TeV 
emission can also be attributed to inverse- Compton scat- 
tering by the same electron population that produces the 
X-ray synchrotron emission. 

Vie wed in X-rays, the young {^^ 330 yr; iFesen et"al] 
|2006[) Galactic SNR Cassiopeia A (Cas A) consists of a 
shell whose emission is d ominated by emis s ion lines from 
O, S i, S, and Fe (e.g. iVink et all [ 1996': 'H ughes et al 
2000[ iWillingale etall |2002l 120031: Hwang fc Laming 
20031 iLaming fc Hwangj 120031 ). Exterior to this shell 



are faint X-ray filaments which mark the location of 
the forward shock. The emission found here is non- 
thermal X-ray s ynchrotron emission from shock accel- 
erated electrons (I Allen et al.l [19971 iGotthelf et all 120011 : 
IVink fc LailliMi 120031 ). " These forward shock filaments 
are observed to exp and with a velocity of ~ 5000 km s~^ 
(jPeLanev fc Rudni ck 2003; Patnaudc & Fcscn 2009), as- 
suming a SNR distance of 3.4 kpc (|Reed et al.lll995() . 

Nonthermal emission filaments are also observed in the 
interior of the SNR and are believed t o be either for- 
ward shock filaments seen in p rojection ([DeLanev et al.l 
120041: iPatnaude fc FesenI |2009| ) or associated with effi- 
cient acceleration of electrons at the SNR reverse shoc k 
(jUchiyama & Aharonian 2008; HelderXvES l2008[ ). 
Fluctuations in both exterior and interior nonthermal fil- 
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aments have also been re ported (IPatnaude fc Fesenll2007l : 
lUchivama fc Aharonianll2008t [Patnaude fc FesenI l2009( ) . 
and the variability is cited as evidence for rapid syn- 
chrotron cooling of TeV electrons in mC-scale fields. A 
two to four year timescale for variations is evidence for 
efficient diffusive shock acceleration in SNR shocks, or 
alternatively the variations are seen as evidence for mag- 
netic f ield fluctuations du e to plasma waves behind the 
shock (jBvkov et al.|[2008l) . 

Emission from Cas A has been seen at en ergies up to ^ 
40 ke V with the Suzaku HXD PIN detector ([Maeda et al.l 
[20091). UP to 100 keV with CGRO OSS E and Integral IBIS 
(|The et al.l[l996l : iRenaud et aT1[2006[ ) , and GeV emission 
has been detected using Fermi-LKT ()Abdo et al.ll2010( ). 
The Fermi observations do not rule out either a leptonic 
origin to the GeV emission from a combination of non- 
thermal bremsstrahlung and inverse-Compton emission 
or a hadronic origin from neutral pion decay. Finally, Cas 
A has been detected at even hi gher TeV energies with 
HEGRA, MAGIC and Veritas ((Aharonian et all I2OOII: 
lAlbert et aO 120071: lHumenskvl l2008[) . Interestingly, the 
centroids for the GeV-TeV emission are located in the 
western region of Ca s A, where the nonthermal X-ray 
emission is brightest (jHelder fc VinkI 120081: iMaeda et ail 
[2OO9I) . 

Here we present Chandra ACIS-S3 observations of Cas 
A taken in 2009 and 2010 which, when compared to 
ACIS-S3 observations taken between 2000 and 2007, 
show the remnant's nonthermal X-ray emission in the 
4.2-6.0 keV band to have decreased at a rate of ~ 1.5- 
2.0% yr~^. In § 2 we discuss our observations, data re- 
duction, and spectral analysis and in § 3 we discuss our 
results and offer some conclusions about the current and 
future evolution of the nonthermal emission in Cas A. 

2. OBSERVATIONS, DATA REDUCTION, AND ANALYSIS 

Cas A has been observed extensively with Chandra 
and we have made use of several GO observations taken 
between 2000 and 2010, including the 2004 VLP (PI: 
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Fig. 1. — Exposure corrected image of Cas A between 4.2-6.0 
keV. The white boxes mark approximately the locations of the 
spectra extracted from the forward shock, while the green boxes 
mark those regions where spectra were extracted from the interior 
filaments. The cyan cross marks the location and 68% confidence 
limit of the Fermi centroid. 



Hwang). We reprocessed each epoch dataset hsted in 
Table [II using CIAC0 version 4.2 and CalDB 4.2.2. The 
more recent observations taken in 2007, 2009, and 2010 
were split due to spacecraft thermal constraints. We 
merged these split observations into a single event lists 
for each epoch. 

We filtered the event lists for those events with ener- 
gies between 4.2-6.0 keV and performed an exposure cor- 
rection assuming a monochromatic 5.1 keV source. The 
exposure corrected image is in units of photons cm~^ s~^ 
pix~^, so to compute the flux of 5.1 kcV photons cm^^ 
s~^, we integrate over the number of pixels. The latest 
dataset, taken in November 2010, is shown in Figure [1] 
As seen in the black curve in Figure [5J there is a clear 
decrease in the overall 4.2-6.0 keV emission from Cas A. 
This decline does not appear to be an instrumental arti- 
fact. Several sources have been observed as Chandra cal- 
ibration targets, including the galaxy cluster Abell 1795. 
In Figure [2l we also plot the 4.2-6.0 keV emission from 
Abell 1795, showing that the emission from this clus- 
ter has not varied by more than 1-2% over ^ 10 years 
(Vikhlinin 2010, private communication). As a further 
check on whether our result for Cas A could be due to an 
instrumental or calibration artifact, we plot in Figure [5] 
the 1.5-3.0 keV emission (exposure corrected at 1.85 
keV) and find that emission in that band has declined by 
less than 1% in 11 years. In addition, the Galactic SNR, 
G21.5-0.9, has been observed extensively with Chandra 
ACIS-S3 as a steady 2.0-8.0 keV continuum source, and 
changes in the flux from that so urce are not observed a t 

e pc _ 

Finally, we note that iKatsuda et al.l ( 2O1O0 reported no 
cha nge in the X-ray sync hrotron emission from SN 1006, 



the level which we re port here (ITsuiimoto et al.ll2010l ). 
ite that iKal 

X-ray sync l 

and iHeinke fc Hoi (|2010l ) reported a decline in the tem- 
perature of the central point source in Cas A, and showed 
that the ACIS-S3 efficiency in the hard band is stable. 

As an additional test of this flux decline, we fitted the 
4.2-6.0 keV continuum emission in Cas A at each epoch 
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Fig. 2. — Comparison of 4.2-6.0 keV flux in Cas A compared to 
the year 2000 observations. The black curve and data correspond to 
changes in the whole SNR, while the red curve and data correspond 
to changes in the western portion of Cas A only. For reference we 
also show the 1.5-3.0 keV flux from Cas A (fluxed at 1.85 keV) as 
well as the 4.2-6.0 keV emission from the cluster Abell 1795. The 
observed decline in the 4.2-6.0 keV emission in Cas A corresponds 
to a fractional decline of -(1.5±0.17)% yr"-"^ across the whole SNR, 
and -(1.9±0.10)% yr~^ in the western limb. 

to a powerlaw model in XSPEC0 version 12.6. We find 
that the modeled flux does decrease with time in this 
energy band, consistent with our analysis of the exposure 
corrected images. As shown in Column 2 of Table [1] the 
fitted spectrum appears to steepen with time. 

The results of our spectral fits listed in Column 2 of 
Table [1] are across the entire SNR and thus can include 
emission from electrons accelerated at the remnant's for- 
ward shock, the reverse shock, or at the contact interface. 
To determine whether the changing emission is from the 
outer forward shock filaments, interior regions or possibly 
both, we extracted spectra from several regions marked 
with boxes in Figure [T] in all datasets and again fit the 
nonthermal emission to a powerlaw model. Galactic ab- 
sorption has a negligible effect on 5.0 keV photons, so we 
did not model it here. 

In Columns 4-6 of Table [l] we list the results from 
the spectral fits for the forward shock regions (the white 
boxes in Figure [1} , the interior regions (green boxes in 
Figured]), and then also just th e interior regions that 
coincide with the Fermi centroid (|Abdo et al.ll201(3 ). As 
can be seen in Table [l] the spectral shape of the forward 
shock filaments varies with time, with Fpg ~ 2.3-2.6 be- 
tween 2000.1 and 2010.8. The interior filaments and the 
region associated only with the Fermi centroid in Fig- 
ure [T] also show an increase in the powerlaw index over 
the same time period. Interior filaments show a steeper 
spectral shape than that seen for forward shock regions. 
We also examined changes in intensity in the western 
portion of Cas A and found that the emission from this 
region is also decreasing with time. This is shown as 
the red curve in Figure [2] and is seen to be steeper than 
the remnant's overall nonthermal emission decline rate 
(black curve). 
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3. DISCUSSION 

As shown in Figure [5] and hsted in Table [U we find 
a decline in the 4.2-6.0 keV emission in Cas A over an 
approximately 11 year span. A linear least squares fit 
to the fractional change in 4.2-6.0 keV emission indi- 
cates a fairly steady decline in the nonthermal emission 
of (1.5±0.17)% yr"-*^ across the entire SNR. In western 
regions coincident with the Fermi centroid, the rate of 
decline appears to be slightly higher at (1.9±0.1)% yr~^. 

A much slower decline in Cas A's radio em ission from 
GeV electrons has been known for some time ( Shklovskil 
[19601: iBaars eFallflQTTl) . iReichart fc Stephens! H2000f ) re- 
port that Cas A has generally been fading at a rate of 
~ 0.6-0.7% yr~^ over a range of radio frequencies. The 
observed decline in Cas A's radio emission has been in- 
terpreted to be a result of the adiabatic expansion of 
the supernova remnant since radio synchrotro n emitting 
electrons have no apprec iable radiative losses ([Shklovskiil 
[19601 lAnderson fc Rudn ick 1996). 

In contrast, the remnant's X-ray synchrotron emis- 
sion is the result of electrons accelerated within the last 
decade, as their radiative loss times are of the order ten 
years, or possibly less. Thus, X-ray synchrotron emis- 
sion is much more sensitive to the present day acceler- 
ation time and radiative losses. For typical SNR pa- 
rameters, synchrotron X-rays are produced in large part 
by the exponential tail of the e lectron distribution (e.g. 
lEllison &: Cassam-Chena^[2005l ). Thus, any energy loss 
will result in a large drop in the emissivity, such as ob- 
served here in Cas A. 

In order to interpret the 4.2-6.0 keV flux changes seen 
in Cas A's X-ray continuum emission, we now examine 
how much of this change can b e attributed to a change in 
the sp ectral properties alone. iZirakashvili fc Aharoniaiil 
(|2007f ) considered the nonthermal X-ray spectrum 
when accounting for losses from diffusive shock ac- 
celeration and synchrotron radiation. In that case, 
the synchrotron emission can be approximated by 
(jZirakashvili fc Aharoniar 



(1) 



N{E) = (t>B'° X E-'°cxp - 




where E is the energy of the emitted photon, which here 
we assume to be 5.1 keV, Ec is the cutoff energy in keV, 
and Fq is the slope of the powerlaw well below the cutoff 
energy. We have split here the normalization constant 
into two factors, for reasons that we will explain below. 
One factor directly corres ponds to the downstreain mag- 
netic field strength (c.f. IGinzburg fc Svrovatskiil I1965L 
equation 3.31), where (j> relates to the overall electron 
acceleration efficiency. Equation [1] implies that the spec- 
tral slope F in X-rays is given by ( Vink et al.,,1999 '): 



F = Fn 



1 E 



2y E. 



(2) 



Thus, the measured spectral index F can be related to 
the cutoff energy by: 

E 



4(F - Fo)2 • 

For Fg, we assume the value inferred from the radio 



spect ral index, namely Fq = a + 1 = 1.78 (|Baars et al.l 
Il977f ). or if nonlinear effects are important, Fq may be 
as low as 1 .25, corresponding to a p article spectral index 
of 1.5 (M alkov fc O'C Drur^^ [20011 ) . For Fq = 1.5, Ec= 
1.0 - 2.0 keV, assuming F = 2.6 - 2.3 and E ^ 5.1 keV. 
Differentiating Equation[2]with respect to time under the 
assumption that only the cutoff energy varies leads to: 



dE, 



Ec dV 
dt \ E dt 



(4) 



According to Table [T] we measure dF/di = 0.018 yr ^ 
in the western portion of Cas A, which translates to 
dEjdt = -(0.032±0.008 - 0.090±0.015) keV yi"^ for 
a.Ec = 1.0 - 2.0 keV. Likewise, from Table[I]for the for- 
ward shock, dV/dt = 0.022 yr~^, which implies dEc/dt 
= -(0.039±0.004 - 0.101±0.008) keV yr'^. 

We can relate the cutoff energy to the shock velocity 
by (jZirakashvih fc AharonianI 120071) : 



Ec w 2.277' 



3000 km s- 



keV, 



(5) 



where r/ is the diffusion parameter, and equals one in the 
case of Bohm-diffusion. Therefore, the change in cutoff 
energy relates directly to a change in shock velocity as: 



dt 



2.0 X 10" 



km s 



The forward shock of C as A has been measured 
to be - 5000 k m s'^ (iDeLanev fc Rudnickl [200l 
iPatnaude fc Fesen|[2009!) . so dVs/dt = -(16±3 - 40±5) 
km s~^ yr~^. The n ature of the reverse sh ock is, how- 
ever, more uncertain. iHelder fc VinS (|2008t ) have argued 
that the western region of Cas A has a velocity as high 
as 6000 km s"'^, implying that the reverse shock in our 
frame is almost stationary. This could indeed be the 
case given a possible interaction w ith a molecular cloud 
in the west ((Keohane et al.l[i996[ ). However, given the 
uncertainties in the nature of the interior filaments and 
the shock velocity of the reverse shock we defer a discus- 
sion on the changes in the western region of Cas A until 
more observations are completed. 

So far we have only interpreted the change in spectral 
slope as due to a change in shock velocity. In addition, 
using the same framework, we can also interpret the ob- 
served flux changes. Differentiating Equation [1] with re- 
spect to time we obtain: 



dN{E) 
dt 



dt 



E- 



exp 



(7) 



where we made use of Equation [4| to rewrite the last 
term. In what follows we will estimate the influence of 
each term separately. 

To start with the last term describing the flux change 
due to the evolution of cut-off energy, it is clear from 
Equation [7| that a change in Ec alone will lead to a frac- 
tional change in flux of: 



1 dN{E) _ dT 
N{E) dt ~^ ~dt' 



(8) 
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This directly connects the change in spectral slope to 
the related change in flux. For the values listed above, 
for both the forward shock and the western region, we 
expect a fractional change of w 4% yr^-'^, which is about 
twice the value that is observed. 

For the other two terms in Equation [7] we must make 
some additional assumptions about their dependence on 
time. Although this adds some uncertainty, it is nev- 
ertheless worthwhile for making order of magnitude es- 
timates. We limit our discussion to the forward shock 
region as both its shock evolution and density structure 
are better understood. 

For the normalization factor 0, we assume that to first 
order it corresponds directly to the amount of gas enter- 
ing the forward shock per unit of time: 



cx ^T^r^pVs (X 



(9) 



where we have used both the fact that the forward shock 
is likely to be moving through the progenitor wind of the 
supernova, hence p cx r~^, and also the idea that the 
forward shock evolution is self-similar with r (x t"^ . Note 
also that this equation is only valid if synchrotron losses 
are large, otherwise (p would be proportional to the to- 
tal amount of material swept up by the forward shock, 
rather than by what is being swept up now. Differentiat- 
ing the above equation with respect to time shows that 
the fractional change in the flux due to a change in is 
estimated to be: 



1 d0 



(m - i)r 



(10) 



For m = 0.66 this corresponds to a fractional flux change 
of —0.1% yr~^, much smaller than the flux changes due 
to a change in cut-off energy. 

For estimating the effect of the second term in Equa- 
tion [7] we make the assumpt ion that cx p Vf , as sug- 
geste d by (iVolk et al.l 120051) . wit h = 2 (jVolk et al.l 
l2005h or B = 3 (IVinljl200l lBeiil[200l . We therefore 
have: 

B^o oc {pVff°/^ cx (i-2™i(™^i)0)ro/2 

for which we have again made use of the idea that the 
density structure falls off as 1 /r^ and r cx . 

The fractional change in B^° , and therefore, the ex- 
pected fractional flux change is expected to be 



1 dB^° 
B^ dt 



^ ^^^-2my.(m-l)/3^)^o/2 



^ ^ro(m(/3-2)-/3)/2 ^ 



(12) 



corresponding to a flux change of -0.45% yr ^ for /? = 2, 
m = 0.66 and Fq = 1.5 or -0.54% yr"! for /3 = 3. These 
values are larger than those due to changes in the first 
term (0), but still smaller than the fractional changes 
predicted from changes in the cut-off energy alone. It is 
important to note that these values are only valid for the 
forward shock, for which we have a reasonable estimate 
of both the preshock density and m. 



There are a number of caveats in the simple analy- 
sis above wh i ch m ay affect our result. For instance, 
iSchure et al.l ()2010D point out that if non-steady state 
acceleration effects are considered, the exponential cut- 
off may have a dependence that is different than that in 
Eq. [2] Additionally, we are averaging over an ensemble 
of shock conditions over a span of ten years; the shock 
conditions and powerlaw spectra will undoubtedly vary 
from position to position, and a sum of powerlaw spec- 
tra do es not yi eld a powerlaw (see also the discussion 
by He lder fc Vin k 2008). In that respect, it is surprising 
that the hard X -ray spectrum of Ca s A is best described 
by a powerlaw (|Renaud et al.l [20 061 ) with little evidence 
for a gradual steepening of the hard X-ray spectrum with 
energy as predicted by Equation [TJ 

It is also possible that not all t he emission aris e s from 
synchrotron emitting electrons. iHelder fc VinkI ()2008l ) 
point out that > 50% of the emission in the west is from 
synchrotron emission, and the remainder is likely from 
thermal continuum. The loss time for thermal electrons 
is much greater than the 2-4 yr timescale seen in the 
TeV electrons. That is to say, over the timescales we are 
investigating, variations in thermal continuum emission 
will be much smaller than variations in the nonthermal 
emission and it is possible that the emission from the 
thermal component may increase over the approximately 
11 yr time span of our observations. 

4. SUMMARY AND CONCLUSION 

Our analysis of the 4.2-6 keV flux of Cas A shows 
a decline of 1.5% yr~^ in the nonthermal X-ray emis- 
sion across the entire SNR over 11 years, with a slightly 
larger decline rate of 1.9% yr~^ from regions along the 
remnant's western limb. We find that qualitatively, the 
observed spectral steepening and decline in fiux can be 
explained by a simple model for changes in the electron 
cutoff energy which are brought about by a natural de- 
celeration of the shock. We estimate an average decel- 
eration of Cas A's forward shock velocity « 10 - 40 km 

The predicted decline in the nonthermal X-ray emis- 
sion is about 4% yr"-'^, which is nearly twice that ob- 
served. The difference between the predicted and ob- 
served decline might be explained by the fact that the 
4.2-6.0 keV continuum emission is not entirely due to 
synchrotron emission from shock accelerated electrons, 
but some of it is from thermal continuum which does not 
evolve on the same timescale as the nonthermal emission. 
We have compared our results to models where the de- 
cline is a natural consequence of either a decrease in the 
number of particles entering the shock or a decrease in 
the efficiency of the shock to amplify the magnetic field, 
and find that these models predict a decline of ^ 0.1- 
0.5% yr~^, which is significantly less than the observed 
decline of 1.5-1.9% yr^^. 

We would like to thank Alexey Vikhlinin and Paul Plu- 
cinsky for useful discussions regarding the stability of the 
Chandra ACIS-S3 quantum efficiency. D. J. P. acknowl- 
edges support from the Chandra GO program through 
grant GO0-11094X as well as support from NASA con- 
tract NAS8-03060. J. V. is supported by a Vidi grant 
from the Netherlands Scientific Organization (NWO). 
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TABLE 1 

Chandra Observations of Cas A and Spectral Fitting Results 



Epoch 

yr 


TsNR 


pa 

10^^" erg cm~^ s"-"^ 




Tlntcrior 


Twcst 


F^'Wcst 

10^^'' erg cm~^ s^^ 


2000.1 


2.81±0.03 


1.61±0.01 


2.32±0.11 


2.66±0.07 


2.66±0.06 


0.229±0.001 


2002.1 


3.01±0.03 


1.56±0.01 


2.43±0.11 


2.75±0.07 


2.74±0.06 


0.223±0.001 


2004.1 


2.99±0.03 


1.54±0.01 


2.42±0.11 


2.70±0.07 


2.73±0.06 


0.215±0.001 


2007.9 


2.98±0.04 


1.45±0.02 


2.55±0.15 


2.70±0.10 


2.80±0.09 


0.197±0.002 


2009.8 


2.99±0.05 


1.42±0.04 


2.61±0.15 


2.78±0.11 


2.78±0.09 


0.195±0.004 


2010.8 


3.07±0.04 


1.34±0.01 


2.56±0.14 


2.82±0.13 


2.85±0.08 


0.183±0.002 



4.2 - 6.0 keV flux 
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